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Abstract. In this paper we study analytically a one-dimensional model for a semiconductor-metal junction. 
We study the formation of Tamm states and how they evolve when the semi-infinite semiconductor and 
metal are coupled together. The non-linear current, as a function of the bias voltage, is studied using the 
non-equilibrium Green's function method and the density matrix of the interface is given. The electronic 
occupation of the sites defining the interface has strong non-linearities as function of the bias voltage due 
to strong resonances present in the Green's functions of the junction sites. The surface Green's function 
is computed analytically by solving a quadratic matrix equation, which does not require adding a small 
imaginary constant to the energy. The wave function for the surface states is given. 
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1 Introduction [T]. Strongly localized surface states are known since a 

long time to be present at the interface of semiconduc- 



The electronic properties of surfaces and interfaces has 
many fascinating features, associated with the formation 
of surface states, modification of the band structure and 
corresponding density of states, and electronic transport 



tors and metals [5] and have been shown to influence the 
conductance of scanning tunneling microscopy (a quasi- 
one-dimensional process) |3]. Also these states have re- 
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cently been found in one-dimensional metamagnetic ma- can be computed analytically is that of a semi-infinite one- 

terials [1]. Therefore the characterization of these states is dimensioiral chain [13J (the surface Green's functions of a 

an important aspect of the physics of surfaces of materials cubic lattice reduce to the solutioir of a onedimensional 

and interfaces between two different solids. problem with analytical solution [22J). 

In the fields of one-dimeirsional physics, specially those Usually, the surface Green's function of a given lead 

related to device applications, the transport properties is computed using a recursive method developed long ago 

are of crucial importante, aird therefore it is important to [14J. The convergence of this method is constrainted by 

study the effect that surfaces and interfaces have on these the value of a small imaginary positive irumber added to 

properties. Of particular interest to us is the non-linear the energy of the Green's function. There is however an al- 

transport through a semiconductor-metal interface. The ternative method to compute the surface Green's functioir 

non-linear transport requires the use of non-equilibrium that does not depend on the value of that small imaginary 

methods, which were first introduced iir this context by part; it even works with the small imaginary part equal 

Caroli et al. [S]. Followiirg these early developments, the to zero, a limit that is implicit in the definition of the 

calculation of the current through a two-band system, con- Green's function jl5[|16| . As we show bellow, this method 

irected to one-dimensional metals, including the effect of can be used to obtain analytical expressions for the sur- 

disorder on the semicoirductor was soon performed [6J, and face Green's function of a system with two atoms per unit 

latter revisited by other authors [7J . As old as these inves- cell, by solving a quadratic matrix equation. 

tigations may be, transport across one-dimeirsional metal- ■ • • ^ ^ j i i- n 

o ^ <^ It IS our purpose, m this paper, to study analytically, 

semiconductor- metal systems is still an active research . . , ^ ii^ir ^- rrn 

withm a simple toy model, the formation ot iamm states 

topic [5], specially in the field of quasi one-dimensional ,,i-,r r -j^ xi-x- j 

f J ^Yie interface ot a semiconductor-metal junction and 

organic conductors. ., ,. , , -c j-- r -n, u 

° its non-lmear transport, as a lunction ot the bias voltage 

Except for very simple models, such as transport across across such interface. We certainly recognize our study to 

an impurity [51 fTUlllllll2| . the transport calculations using be that of a toy model problem, but the fact that all quan- 

the non-equilibrium Green's functions method are all nu- titles can be computed analytically makes our results rele- 

merical. This is so due to the fact that both the surface vant for the study of more sophisticated models where no 

Green's function of the contacts and the Green's func- analytical solution exists. We will compute all the needed 

tions of the system (often called device) are obtained by Green's function analytically, which in turn give analyti- 

the inversion of very large matrices, a process that in gen- cal expressions for the tunneling probability and for the 

eral has no analytical solution. There is even the wide- system density matrix. To our best knowledge this study 

spread idea that the only surface Green's function that has not been done before. 
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The paper is organized as follows, in Sec. [5] we formu- 
late the problem in terms of an infinite dimensional block 
tridiagonal matrix and explain how to derive the needed 
Green's function as well as the method of solution for the 
surface Green's functions. In Sec. [3] we study the existence 
of surface states in the semiconductor and metal surfaces 
and how these states evolve to sharp resonances when the 
semiconductor and the metal are coupled. In Sec. 2] we 
compute the charge current through the interface for arbi- 
trarily large bias voltage and the local electron occupation 
at the sites of the interface. Finally, in Sec. [S]we give our 
conclusions. 

2 Hamiltonian toy model and formalism 

In order to give a analytical description of the properties 
of a semiconductor-metal junction we consider here a one- 
dimensional model. In Fig. [T] we depict the model used. 
The unit cells have two atoms. In the case of the semicon- 
ductor the atoms are different (a and b in Fig. [1]); in the 
metal all the atoms are equal (c in Fig.[T]). The doubling of 
the unit cell in the metal is only a matter of convenience, 
since it allow us to write the Hamiltonian as a tridiag- 
onal block matrix (see below). The junction takes place 
in the unit cell n = 0. The terminal atoms of the semi- 
conductor and the metal are represented by the letters 
a' and c', respectively. In the semiconductor the hopping 
within a unit cell is J and between unit cells is V. In the 
metal the hopping is — t. The local energies in the semi- 
conductor are and for atoms a and 5, respectively, 
and for the metal atoms the local energy is represented by 
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Ec- In geireral we expect that the atoms at the surface of a 
metal, say, have different local energies and be comrected 
to the bulk atoms by different hopping parameters. This 
fact introduces the parameters e^, e^, V' , — t'. The hop- 
ping across the interface is W. The reader may note that 
the model has a somewhat large irumber of parameters. 
This is necessary, since at the interface of both the metal 
and the semiconductor both the hopping and the on-site 
eirergies are modified relatively to their bulk values. Nat- 
urally, the number of parameters can be greatly reduced if 
we do not take into consideration the effect of the surface 
on the model parameters. Although our study is aimed 
to a general description of the effects of surfaces on the 
non-linear transport across an interface, its application to 
a real system is a possibility; this would require us to use 
the values of the parameters appropriate to the system 
under study. We note, however, that all energies are given 
in units of t, and therefore the parameter values used in 
the simulations can apply to a general system. 

J V J V J V w -f -t -t 

^^^^-<6^ 0—0—0—0— 

-3 -2 ""o 7" 

Fig. 1. (color online) Model of a semiconductor-metal inter- 
face. The parameters we used through the text are: Ec = 0.1, 
J = 0.8, V = 0.7, ea = 0.1, et, = 0.2, t' = 0.9, V = 0.9V, 
= l.lea, Ec ~ S.Oec, where all energies are in units of t. The 
parameter W will be varied, but is also given in units of t. 

If we denote the semiconductor Hamiltonian by H^, 
the metal Hamiltonian by Hj^, and the junction Hamilto- 
irian by Hd (defined by unit cell n = 0), the full Hamil- 
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tonian of the problem, written in the local atomic basis, of the system is defined by 
reads 



H 



Hl Vl 

yl Hd vl 

VrHr 



{lE + iQ+ -H)G+ = 1, 



(7) 



(1) where we have chosen the retarded function (denoted with 
the -|- superscript), and 1 is an infinite identity matrix. 



where the matrices Vl and Vr are the coupling between ^^^^ ^j^^ Green's function is 

the semiconductor and the junction and the metal and 
junction, respectively. If we write and Hr explicitly, 
the Hamiltonian ([T]) acquires a block tridiagonal form and 
reads 



H 



C Hgm Al 



A\ Hd Ar 



/it, Hm Cr 



C\, H„i Cr 



Gll Glci Glr 

GdL Gdd GdR ■ (^) 

GrL GRd Grr 

The quantity of interest is Gddi which is shown to have 
the form 



Gld = {^E + iQ+-Hd-St-S+)-\ 



(9) 



(2) 



where the matrices aird are the self energies and 
have the form 



I^t = A{GlLAL, S+ = ArG+^A\ 



R ' 



(10) 



where the several matrices are given by 



Ham. — 



Hm — 



Ar^ 



and 



fa J 
J £6 

-t Cc 


V 



Hd^ 



Gr = 



Ar^ 








4 w 
W 


V 



-t' 



(3) 



where the Green's functions Gj^ aird GJ^^^ are the sur- 
face Green's function of the Hamiltonians Hl and Hr^ 
respectively. These Green's functions are defined as 



(4) 



Gr 



Gj^ Usm Gl 



G\^ U sm 



(11) 



• • G-2,-2 G-2,L 

■ ■ Gl,-2 Gll 
with Usm = E+iO^ — Hsm and a similar equation defining 
(5) G^^. It is possible to find a close form for Gj^ and G 
P^lfT^. reading 



RR 



(6) 



G+^ = [E + zO+ - Hsr,, - gIGllGl]-' , (12) 
G+^ = [E + iO+ -Hm- GrGrrG^]-' . (13) 



Since we want to study the properties of the junction, 
we ireed to compute local quantities. This is best accom- (Iir the Appendix we give a simple derivation of Eqs. (|12p 
plished using Green's functions. The full Green's function and ([12]).) The solution of and can in general be 
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done numerically only, by using a decimation procedure Equation (|16p has non-trivial solutions if the following de- 
[13] , or a direct iterative solution [T3] . In these two meth- terminant is zero 

ods it is necessary to introduce a small imaginary part, ,,^+9 ^ „j_n ^ ,, „ 

l|Cl?/,^ + (ffsc--B-20+)y, + Cl|| =0. (17) 

that is 0^ is replaced by ?7+, where is a finite number. 

The solutioir of Eq. (|17p produces in principle four eigen- 
values yi, of which only two are physical. In general the 
solution of the quadratic matrix equation has to done nu- 
merically, but for our 2x2 matrix an analytical solution 
exists. For real eigenvalues yi. convergence of the solution 
for G\i^ requires that iji <\. All the imaginary eigenval- 
ues satisfy the condition \yi\ = 1, that is the imaginary 
yi's can be written as 



The rate of convergence of the two methods depend of 
the value of 77, which we would like to be as small as pos- 
sible. There is however another method available which 
is based on the solution of a quadratic matrix equation 
P^[TC1[T7| and that does not require the use of a finite 
value of r\. This method was recently used in the context 
of transport through molecular junctions |18II19| . but has 
been essentially forgotten by the community working on 
surface Green's functions applied to non-linear transport. j/i = x ± — , x <1, (18) 



Let us explain how this last method works considering, for ^nd the choice of the sign is made in order to satisfy the 
this purpose, the solution of G\^. We start by defining an analytical properties of G+^, namely 5G+^ < 0. Aher- 
auxiliary quairtity Y = G+^Cl, which allows to write Eq. natively, a small positive imaginary part can be added to 



m as 



gIy^ -f {Hsc -E- iO+)Y + Cl=0. (14) 



the eirergy, and the correct choice of yi are those solutions 
that lie in the unit circle. 



We now assume that there is a similarity transformation oi_ _ ^- ^ -i-u- 

^ 3 Local electronic properties at equilibrium 

Q that diagoiralizes the matrix Y, defined as y = QYQ^^, 



where F is a diagonal 2x2 matrix, Y — diag(?/i, 2/2)- Using 
Y we can write Eq. (fT4|) as 



In this section we want to address the electronic properties 
of the junctioir when the chemical potential of the semi- 
coirductor and the metal are equal, and therefore there is 
Gf^QY + {Hsc — E — iO^)QY + GlQ — . (15) current flowing through the system. To that end we 

need to compute Gj^, G^^, and G^^; they are all 2x2 
matrices. For the case of Gj^, Eq. (fT7)l has the form 



Writing Q as Q ~ ((71,(72)7 where qi and (72 are vector 
columns of two elemeirts, we obtain that the solution for 
Y and Q reduces to the solution of a quadratic eigeirvalue ViiE - <^a){E - e^) - yiiJyi + V){J + Vyi) = , (19) 



problem of the form which has three solutions, yi ~ yi — and 

[Glyf + {Hsc-E-iO+)y, + GL]q^^O, i-l,2.(16) y, = ^ ± y/ ~ 1 , (20) 
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with p =[{E- ea){E - Eb) - J2 - V^]/{2VJ). For /? > 1 with a ^ [{E - e^f - 2t^]/{2t^). 



the correct choice of yi is 



y.-Z/2=/3-sgn(/?)v//32-l. 
and for /3 < 1, the correct choice is 



(21) 



The calculation of Gj^ requires the determination of 
the self energies. These are simply obtained as 



r+ — 



(30) 



For 2/1 the eigenvector q\ is 





(22) and 



qi 



Vl 



(31) 



(23) 




.Ot''(G+;v)ii_ 

The matrix elements (6*^^)22 and (G^^)ii are, after some 
algebra, simply given by 



where vi in any real number. For y2, the eigenvector 52 is 

(24) 



Q2 



U2 
V2 



(0+^)22 = iE-et,)-\l + Jy2/V). 
iG+j,U = {E-e,)-\l + y2). 



(32) 
(33) 



with U2 — V2X, and X given by 

J + y2V y2{E-tb) 



X = 



(25) 



At first sight, Eq. does not look like the surface Green's 
function of a semi-infinite one-dimensional chain (that is 
because we used two atoms per unit cell for the metal) 



E-ea Jy2 + V ' 

with V2 any real number. It is now a simple task to com- f^' however simple algebraic manipulations show that 



pute Q and its inverse, from which Y, and therefore Gj^, 
is obtained. The surface Green's function is given by 

E-€b J 

J E-(a- Vy2/X 
Using exactly the same procedure we obtain for Gjjj the 
equation 



(Gj^jj)ii can be put in the known form 

E-Er 



(34) 



y2 

JV 



2t^ 2^2 

(26) foi' {E — Ec)^ < and a similar equation for (E — > 



4^2. Using Eqs. ^ and Gj^ is given by 



G 



dd 



S{E) 



G 



RR 



t2 



with 



and 



Z = 



E-ec + y2tZ -t 
-t E-e, 



E 



SciE) W 

W Sa{E) 



(35) 



(27) with S{E) ^ SaiE)SciE) ~ W^, and 



til+y2) ^ 

y2{ec-E) t(l-fj/2) 



2/2 



a — sgn{a)yc?~l , > 1 , 

a — i\J\ — oP- , < 1 , 



(28) 



(29) 



Sa{E)=E-e[,-V'^{G+^)2 
S,iE)^E-e',-t'\G+j,)n 



(36) 
(37) 



If = 0, the two systems are decoupled, and (Gj^)ii 
is the surface Green's function of the semiconductor and 
(^(L)22 is the surface Green's function of the metal. 
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In some conditions, the existence of surfaces in a ma- 
terial give rise to surface states. In semiconductors these 
states lie in the gap of the semiconductor. These type of 
states are determined by the condition S{E) ~ 0, with 
S{E) a real number. When S{E) is imaginary, S{E) = 
Si{E)+iS2{E), resonant states are determined from Si{E) 
0; 5*2 (-E) will be related to the width of the resonance. The 



arrosfi a sem 



ir.ondiirt or- metal innrtinn 




local density of states at the iunction atoms is given by y.,. „ / , i- \ r>i ^ j-ii • ■ /j i j i- n j 

iig. 2. (color online) Plot or the imaginary (dashed line) and 



(either with zero or finite W) 



Pa' 

Pc' 



TT 



the real parts (solid line) of Sa(E) and SdE). To make the 
resonance in the metal well visible we have used ej. = See- In 
the lower panels we have the local density of states at atoms a' 
(39) and c' when W = 0. The other parameters are those introduced 
in Fig. [U 



(38) 



where pa' and pc' are the local density of states at a' and 
c' atoms, respectively. 

Considering first the case W — 0, the case in which the 
two systems are decoupled, the surface states of the semi- 
conductor satisfy the condition Sa{E) = 0, with Sa{E) 
real, and the resonances satisfy the condition ^Sa{E) = 0. 
Similar expression hold for the metal with Sa{E) replaced 
by Sc{E). In Fig. [5] we plot the imaginary and the real 
parts ofSa{E) and SdE). 

In the case of the metal, the condition diSc{E) — cor- 
responds to the maximum of the density of states, which 
represents a very broad resonance. In the case of the semi- 
conductor we see that in the gap we have Sa{E) = 0, with 
Sa{E) real, thus corresponding to a pole in the Green's 
function and, therefore, represeirting a surface or Tamm 
state. Let us calculate the wave function of the Tamm 
state (for simplicity we consider the case ea' ~ 0). The 



general wave function has the form 
-1 

= ^ {a„i\m,a) +brn\m,b)) . 



(40) 



where |m, a) and |m, b) are position-base states at the sites 
of the chain, and a„i and 6,„ are the corresponding ampli- 
tudes. The matrix elements of the Hairriltonian, (n, a|iJ|?/') 
and {n,b\H\ijj) leads to 

V'b-iS-i,n + Vbn9{-n + Jbn+l = Ettn , (41) 

Ja„_i + Van = Eb„ , (42) 

where 6{x) is the Heavyside function, with 9{0) ~ 0. The 
above equations are subject to the boundary condition 
bo = 0. We now make the observation that for E ~ 0, 
there is a non trivial solution for the amplitudes a™ given 
by the recursive relation 



Jar, 



+ Va^ 



0, 



(43) 
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and all bm — 0. The wave function of the surface state is 
therefore given by 

-1 

|^)= ^ ai(-l^/J)l-l-iKa), (44) 

m— — oo 

with ai = -y/l — (V/J)^ obtained from the normalization 
of the wave function. Clearly, for the surface state to exist 
we need V/J < 1 , the case we used in our numerical calcu- 
lations. For V/ J > 1 the surface state is absent. In the par- 
ticular case V/ J > 1 a Shockley state will develop in the 
gap, which is not a surface state. Also for the metal there 
are some conditions where localized states can form at 
the surface. Let us take the semi-infinite one-dimensional 
metal introduced in Fig. [TJ whose Hamiltonian reads 

oo 

H = -tY,{\n >< n + 1\ + H.c.) + |0 >< 0|ec' , (45) 

and we have assumed all the hoppings equal. Proposing a 
localized wave function of the form 

oo 

^ioc = A^e-"^|n>, (46) 

Tl = 

and writing the Schrodinger equation as 

- tai ^ {E ~ t,,)aQ , (47) 
Ean = -i(a„_i + a„+i), n > , (48) 

we obtain for the energy of the localized state 

E = -2tcosh\, (49) 

with A the solution of 

e^^ec'/t. (50) 

Since A must be larger than zero we must have Cc' /t > 1 
and the energy of the localized state is 

(51) 
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located below the bottom of the metal band. So, in this 
special coirdition it is possible for the metal to develop a 
localized state at the surface. If we generalize the above 
case and include the possibility that the hopping between 
the site c' and the site c is i', the wave function of the 
localized state has to be generalized to 

oo 

V;ioc = A^e-"^|n>-|-A6|0>, (52) 

n=l 

and the Schrodinger equatioir has now the form 

- t'oi = (£■ - ec')ao , (53) 
Eai = -t'ao + tai , (54) 
Ea„ = -<(a„_i + a„+i), n > 1 , (55) 

whose solution gives b = t/t' and 

2e"' = ^#^±^7T^v/(^^^0^+W^, (56) 

and the energy is still giveir by i? = — 2tcoshA. Since 
we must have A > only some values of the parameters 
produce surface states, in particular we must have t' > t 
(which is not the case considered in the simulations.). 

We now make W finite coupling the semiconductor 
and the metal. In Fig. [3] we plot the imaginary aird the 
real parts of S{E) for different values of W. Because S{E) 
contaiirs now contributions from both the imaginary parts 
of Sa{E) and SdE) there will be a finite imaginary part 
in the energy range of the semiconductor's gap. As a con- 
sequence the Tamm state previously located in the gap 
becomes now a sharp resonance aird is visible in the local 
density of states pa' , as we show in Fig. Of particular 
interest is the strong transfer of spectral weight from the 
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Fig. 3. (color online) Plot of the imaginary and real part of 
S{E) for different values of W . The parameters are those in- 

j J J r?:™ Q 




Fig. 4. (color online) Plot of local density of states p^/ and 
Pc' for different values of W . The left panels depict the semi- 
conductor density of states at site a' and the right ones the 
same quantity for the metal at site c'. At site c' the formation 
of an anti-bound state is clear, where a dip in the density of 
states is seen. The parameters are those introduced in Fig. [T] 

density of states of the metal to the resonance in the semi- 
conductor, given rise to an anti-resonance in the density of 
states of the metal. It is interesting that the high density 
of states persists where the resonance was located even 
when ^S{E) is no longer zero, albeit small. Only for large 
values of W, such as = 1, a large transfer of spectral 
weight to resonances at lower and higher energies takes 



place. From Fig. [3] we see that for moderate values of W, 
even when we have ^S{E) = 0, no sharp resonance appear 
in the local density of states because QS{E) is very large 
at those points. For W ~ 1 strong resonances appear at 
the lower and top edges of the bands. 

4 Transport properties 

We now study the non-equilibrium transport across the 
junction. This is done using the non-equilibrium Green's 
function method |21| . This method is particularly suited to 
study the regime where the system has a strong departure 
from equilibrium, such as when the bias potential Vf, is 
large. The system is however in the steady state. Since the 
seminal paper of Caroli et al. on non-equilibrium quantum 
transport [5j, that the method of non-equilibrium Green's 
functions become generalized to the calculation of trans- 
port quantities of nanostructures. There are many places 
where one can find a description of the method [^1231 
I24| . but an elegant one was recently introduced in the con- 
text of transport through systems that have bound states, 
showing that the problem can be reduced to the solution 
of a kind of quantum Langevin equation |25| . 

The general idea of this method is that two perfect 
leads are coupled to our system, which is usually called the 
device. In our case the device is defined by the junction, a 
two site system, involving the a' and c' sites. The Green's 
function of the device has to be computed in the presence 
of the leads. This corresponds to our Gj^ Green's function. 
Besides the Green's function we need the effective coupling 
between the leads and the system (the junction) which are 
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determined in terms of the self-energies as 

(57) 



i rrin H 1 1 nr- m a 1 i 1 1 n r t i on 



W=0.1 


1 ' 1 ' 


-- W=0.25 




-- W=1 


' f 


- W=2 


1 ' \ 




1 1 \ 









therefore the effective coupling /l/h depends on the sur- 
face Green's function of the perfect leads. According to the 
general theory the two leads are in thermal equilibrium at 
temperatures Tj^/fj and chemical potential (J-l/r and are 

connected to the system at some time in- The bottom line _ , , ,. \ m . . mri-,\ ^ tr „ , n-r 

Fig. 5. (color onhne) Transmission T(E) tor \/ = and dit- 

is that the total current through the device is given by ^^^^^^ ^^^^^^ ^ ^^^p j^^^^. ^^^^ ^.^^^^^^ ^^^^^^ 





(both spins included) 



J 



2e 



dST(£;) [/(£;, Mi , Ti) -/(£;, ^i^, T^)] , (58) 
where the transmission T{E) is given by 

T{E) = A^^Tr[rLG+,rRG-J . 
Performing the trace we obtain 



(59) 



and W — 0.25 (down left). Non-linear current J as function of 
the bias energy eV (top right), for different values of the chem- 
ical potential and W — 0.25. Differential conductance (bottom 
right), dJ/d{eV), for the same curves in top right panel; the 
negative values of the differential conductance are due to the 
non-linearities of the current curves upon the bias potential. 
The parameters are those introduced in Fig. \T\ 



T{E) ^ W'H"^{G+^)22'^iG+^U 



\sm' 



(60) 



Because we are applying a bias voltage across the junc- the discrete Poisson equation ignoring the charge fluctua- 
tion, we choose = /i and fiR ^ fJ. ~ eV. Also the elec- tions taking place to screen the external electric field. (In 
trostatic potential has to change continuously between the general terms, to determine the potential at sites a' and 
semiconductor and the metal. We choose that the varia- C we would have to solve the Poisson equation coupled 
tion of the potential is proportional to the distance to the the solution of the Schrodinger equation EHIEZ], but 
electrodes, therefore the local energies in the junction have the above transformation of e^' and e^' gives a good first 



approach to the exact result.) One technical aspect worth 
stressing here is the fact that, in general, non-linear trans- 
port should be done in a self-consistent way. This is the 
The sites of the metal are shift by —eV. Since T{E) is case in two situations: (i) when interactions (Coulomb or 



to be modified according to 

ea' ^ ea' - eV/3 , e,- ^ e^' - 2eV/3 



(61) 



computed using and this depends on Ca' and on ec', phonons, say) are taken into account; (ii) when the poten- 
T{E) will also be a function of the bias potential V. The tial inside the conductor is relevant for transport. In the 
choice made in Eq. (1611) corresponds to the solution of case we are considering here the conductor is reduced to a 
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two-site system in the absence of interactions. Therefore, 
a self-consistent calculation is not needed. 



W=0.25, (1=0.2 



Since the transmission function T{E) depends on the 
imaginary part of the surface density of states, the gap 
due to semicoirductor shows up. From Fig. [5] we see that 
the current J across the junction is only linear upon V 
for very small values of the bias potential. As V is fur- 
ther increased non-linearities in J start to develop due to 
the energy dependence of T{E). which is strongly influ- 
enced by the resonances shifting in energy as we vary the 
bias (see left down panel of Fig. [5]); for large V (say ^ 
2) the upper band of the semiconductor non longer con- 
tributes to the transport, due to the relative shift of the 
local density of states induced by V . Also at large V reso- 
nances develop in T{E) at both lower and higher energies; 
this effect is seen at the bottom left panel of Fig. [S] for 
V — 2.4 eV. As the junction is strongly biased the current 
starts being suppressed due to the relative motion of the 
local energy of sites in the junction and the bands of the 
semiconductor. For large bias, the on-site energy of the 
junction's sites moves to lower energies and these states 
are no longer available to the electrons coming from the 
semi-conductor. A qualitative model for tunneling across 
a one-site system also presents the general features seen 
for J in Fig. [51 except that the effects due to resonances 
are not present [53]. Naturally, when the chemical poten- 
tial of the semiconductor lies on the gap it is necessary a 
finite bias voltage to produce a current. 
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Fig. 6. (color online) Top figure: occupation numbers n^,/ and 
as function of the bias voltage V for different values of 
the chemical potential. Bottom figure: matrix elements of the 
Green's functions: Gil = -yLUG+Jn; G12L = 7l|(G+ )i2l^; 
G12i? = 7fl|(G]'Ji2p; G22R = 'yR\{G+^)22\^ . The parameters 
are those introduced in Fig. [1] 

The average number of electrons, per spin, at the site 
j of the device is given by 

/oo 
dE{G+^r^G;i,),,f{E,pi^,n). (62) 

Specializing to the sites j = a' and j — c' we obtain 

na' = I di?[7L|(G+ )iiP/l + 7fll(G+ )i2p.fe] , (63) 
rf^^[7L|(G+ )i2|Vl +7fll(G+ )22|Vfl.] , (64) 
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where the Jl/r are the Fermi functions, and jl/r a-re proportional to VF^ . So, in this case, the occupancy of the 
defined as sites a' and c' is only due to the electronic wave-function 



y/2 coming from the system to which the corresponding site 

7L - - — S(G+j22 , (65) 

,2 is directly connected. 

IR - -—^{G^hi ■ (66) 

For large bias voltage the matrix elements that con- 
From Figure E we see that the local electronic deir- tribute to n^' develop large resonances which contribute 
sity at a' and c' sites is a highly non-linear function of to the increase of relatively to n^' at F ~ 3 eV. 
the bias potential. This behavior can only be understood It turns out that the details of the behavior of Ua' 
by looking at the behavior of the matrix elements of the and n^' depend some what on the relative strength of the 
junction Green's function both as a function of energy hoppiirg parameters and oir-site energies, 
and bias voltage. Clearly there are many resonances in 

the Green's function matrix elements, which renders the 5 Summary and conclusions 

analysis rather difficult. For small and intermediate bias 

In this paper we have studied a simple one-dimensional 

voltage we have: 

model of semiconductor-metal junction. The advantage of 
'YLliG^cdi-i 7-R 1(^^^)121^ , (67) this simplification is that all the features can be stud- 

ied using analytical expressions. We have shown that the 

and 

semi-infinite metal does generate surface states in partic- 

7L|(G+)i2p»7fl,|(G+)22r (68) 

ular conditions. For the semiconductor surface states can 

These inequalities are due to the relatively large value of r . ti i r j 1 1 j -i 

^ J t> form m the gap. Kesoirairces can be lormed local density 

e' relatively to e^. On the other hand, when the voltage r . . r , 1 ^ 1 -r ii r i 
c J c , *= 01 states of the metal it the surface parameters are very 

increases, the above inequalities transform, roughly, into i-rr , r ,1 • u n rpi •.• r 

' ^ I to J J dmerent irom those m the bulk, ihe energy position 01 

approximated equalities. , ■ ■ u 

^ ^ such resonance is given by 

As a general trend, and rir' will decrease at rel- „ ,„ 

= 0/2 m ■ (69) 

atively large values of eF, due to the dependence of the it — i 

energies ta' and t^' on V (see Eq. ((6T|) ): this is specially In the case t' = i and ec' — Cc, Er is simply the energy of 

the case for small coupling between the semiconductor and the maximum of the local density of states, 

the metal (small W\ In the case of small coupling, the When the interface is formed by making the param- 

contributioii due to the off-diagonal Green's function is eter W finite, the surface state formed in the gap of the 

small, since this latter quantity is proportional to W and semiconductor becomes a resonance, because the density 

therefore its contribution to n^' and ric' (see Eq. (j64p ) is of states of the metal is finite in the gap of the semiconduc- 



The Hamiltoiiian of the semi-infinite lead has the form 
Hr Cr 



H 



Cj^ Hr Cr 

Cr Hr, Cr 



Hr Vr 



N. M. R. Peres: Analytical study of non-linear transport across a semiconductor-metal junction 

tor. This is a consequence of the choice of the parameters A Simple derivation of Eq. ( I13p 

for the semiconductor for the metal. We could as well have 
chosen a different set of parameters, such that the density 
of states of metal was zero in the gap of the semiconductor. 
In this case the surface state, as pole of the full Green's 
function, would still survive as long as the renormalization 
of its energy due to the finiteness of the real part of the 
metal Green's function would not move it away from the 
gap. As W changes the resonances in the semiconductor 
and in the metal drift in energy 

These resonances have a strong impact on the trans- 
port properties of the junction, because for fixed values 
of W, their energy location is dependent on the bias volt- 
age which changes the on-site energies of the atoms in 
the junction. Their effect is even more dramatic on the 
electronic occupancy of the atoms in the junction. 

Although our system is a very simple one, the features 
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(70) 



Let us define the Green's function of the system by {E — 
H)G ~ 1. The Grr Green's function is obtained from 



E-Hr -Vr 



-V^ E-H 



Grr Gri 
GiR Gil 



(71) 



From Eq. ([7T|) we derive 



{E — Hr)Grr — VrGir — 1 , 
-V^Grr + {E- H)GiR = . 

Solving this linear system for Grr we obtain 

[E - Hr)Grr - Vr{E - H)-^vIGrr = 1 , 



(72) 
(73) 



(74) 



seen in this case will also be present on more reahstic but {E - H)-^ = G and VrGVI = CrGrrC'^^, which 
cases. However, in the context of quasi one-dimensional leads to Eq. 



organic conductors our calculations have direct relevance; 
for example, the electronic transport in a molecular wire 
will develop features as those described here close to the 
contacts to the metallic leads. 
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